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ABSTRACT 

We present Atacama Large Millimeter/submillimeter Array (ALMA) observations from the 2014 Long Base¬ 
line Campaign in dust continuum and spectral line emission from the HL Tau region. The continuum images 
at wavelengths of 2.9, 1.3, and 0.87 mm have unprecedented angular resolutions of 0."075 (10 AU) to 0."025 
(3.5 AU), revealing an astonishing level of detail in the circumstellar disk surrounding the young solar ana¬ 
logue HL Tau, with a pattern of bright and dark rings observed at all wavelengths. By fitting ellipses to the 
most distinct rings, we measure precise values for the disk inclination (46.72° ±0.05°) and position angle 
(+138.02° ± 0.07°). We obtain a high-fidelity image of the 1.0 mm spectral index (a), which ranges from 
a rsj 2.0 in the optically-thick central peak and two brightest rings, increasing to 2.3-3.0 in the dark rings. The 
dark rings are not devoid of emission, and we estimate a grain emissivity index of 0.8 for the innermost dark 
ring and lower for subsequent dark rings, consistent with some degree of grain growth and evolution. Addi¬ 
tional clues that the rings arise from planet formation include an increase in their central offsets with radius and 
the presence of numerous orbital resonances. At a resolution of 35 AU, we resolve the molecular component 
of the disk in HCO + (1-0) which exhibits a pattern over LSR velocities from 2-12 km s -1 consistent with Ke- 
plerian motion around a ~1.3M© star, although complicated by absorption at low blue-shifted velocities. We 
also serendipitously detect and resolve the nearby protostars XZ Tau (A/B) and LkH<r358 at 2.9 mm. 

Subject headings: stars: individual (HL Tau, XZ Tau, LkHa358) — protoplanetary disks — stars: formation 
— submillimeter: planetary systems — techniques: interferometric 
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1. INTRODUCTION 


Inside the Taurus star-forming complex, HL Tau is a 
young star located in a molecular ridge of length ~ 0.05 pc, 
which forms part of the wall of a large-scale bubble seen 
in both 13 CO and faint scattered light (IWelch et al.l 120001: 
Angl ada et al.l 120071) . This material contributes to the esti¬ 
mated ^ 24-33 magnitudes of visual extinction (Ay) towards 
HL Tau (IClose et al.lll997fc IMen’shchikov et al.lll999l) . Three 
other less obscured protostars lie toward the edges of this 
molecular ridge, in order of separation from HL Tau: XZ Tau, 
LkHa358, and HH30 (e.g.. iMoriartv-Schieven et al . 2006 ). 
Together these four protostars make up the “HL Tau region” in 
LI551. In this paper we adopt the standard mean distance to 
Taurus of 140 pc (IRebull et al.ll2004l) for the HL Tau regioi0 


As a result of the high extinction, the central star of 
HL Tau has not been directly detected at optical wavelengths, 
and only a conical refl ection nebula has been observed (e.g. 
Stapelfel dt et al.lfl995l) . Near-infrared (NIR) images reveal a 
point-like object that has been attributed to direct stellar radia¬ 
tion; though, scattered emission from a central hot disk is also 
a like ly explanation (IClose et al.l 119971: [Men’shchikov et all 
119991) . Based on high dispersion optical spectroscopy, HL Tau 
is classified as spectral type K5±l ([White & Hillenbrandl 
2004 ). Analyses of the full SED find bolometric luminosities 
ranging from 3.5 to 15 L 0 , and classify it as a Class I-II pro¬ 
tostar with evidence of an extended envelope and a circum- 
stella r disk ([K enyon & Hartmannl ll995l: IMen’shchiko v et all 
1999; Robitai lle et al.l 120071) . In all, the stellar properties of 
HL Tau still have a high degree of uncertainty despite exten¬ 
sive study. 

Measurements in the millimeter regime suffer far less from 
extinction effects, and the HL Tau disk is one of the bright¬ 
est at these wavelengths ([Andrews & Williams!120051) : hence 
HL Tau has been a favored interferometric t a rget over the 
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a few). The highest angular resolution mm observations to 
date have a 0."13 (18 AU) beam at 1.3 mm in which the 
disk is resolved with an outer radius of 120 AU at posi¬ 
tion angle +136° and an inclination (from face-on) of 40° 
dKwon et al. 2011 ). Modeling of previous millimeter data 
suggest an HL Tau d i sk mass of Mj y 0.0 3 -0.14 M^ 
(Rob itaille et al.l 120071: iGuilloteau et al.l 120111: iKwon et ail 
1201 ll) . This high disk mass is within an order of magni¬ 
tude of previous HL Tau stellar mass estimates ranging from 
0.55 to 1.2 M^ ([Sargent & Beckwithll99HIClose etal.il 19971: 

Whit e & Hillenbrandll2004l) . suggesting that the disk may be 

close to being gravitationally unstable. 

At longer (cm) wavelengths, free-free emission from an 
outflow jet ortho gonal to the dis k plane contributes to 
the continuum (IWilner et all 119961: iRodmann et al.1 [2006; 


iGreaves et al.l 120081: ICarrasco-Gonzalez et al.l 120091) . Also 
known as HH150-151, this jet is seen in optical and near- 
IR shock-tracing lines such as [SII] and [Fell], extending 
up to r sj 10,000 AU from the star at vel ocities of 100-200 
kms -1 dMundt et al.ll990tlkrist et al.l2008l) . The inner jet has 
been traced down to ~ 100 AU from th e disk dPvo et al 1 12006: 
iTakami etaH 120071: iBeck et all 120101) . At mm wavelengths, 
CO emission traces a lower-velocity, entrained component of 


20 However, Very Long Baseline Interferometry measurements indicate 
that the line-of-sight depth of the Taurus complex is ~20 pc, so that reported 
linear distances could be in error by up to 15% (e.g., ILoinardll2013l) 


this outflow, mainl y showing red-sh i fted gas on the sou thwest 
side of the disk (ICabrit et al.l 1 1996t iMonin et all 1 19961) . The 
central few hundred AU of this outflow was also observed 
using the SMA in the 12 CO (3-2) line dLumbreras & Zapatal 
120141) . and showed a similar dominance of red-shifted gas. 
Detailed study of the molecular gas associated with the disk 
has heretofore been plagued by spatial and kinematic confu¬ 
sion due to the extended ridge and outflow emission near the 
ambient velocity of vi sr = 6-7 km s -1 (e.g. ICabrit et al J-19961) . 

Overall, HL Tau is an excellent example of a system just 
emerging from its protostellar cocoon (i.e. evolving from 
SED Class I to Class II), which contains a massive com¬ 
pact disk as well as highly collimated outflow. The sys¬ 
tem is youn g (<l -2 Myr based on the cluster age of Tau¬ 
rus, [BriceIo_eral]|20Q2), and the high disk mass makes it an 
ideal candidate for disk insta bility and early planet formation 
(e.g. llNero & Bjorkmanl2009l) . In this paper, we present multi¬ 
wavelength ALMA Science Verification (SV) continuum and 
spectral line observations from the 2014 Long Baseline Cam¬ 
paign (LBC) of the HL Tau region. These data dramatically 
demonstrate the revolutionary impact that the full sensitivity 
and resolution of ALMA will have on the field of star and 
planet formation. In this work, we aim to give an initial taste 
of the incredible richness extant in these first ALMA long 
baseline observations of a protoplanetary disk, and we expect 
that in time the community will uncover the full scope of what 
is possible with these data. 


2. OBSERVATIONS AND DATA REDUCTION 

The ALMA SV data on HL Tau was taken between 2014 
October 14 and November IA with baselines ranging from 
15 m to 15.2 km. lALMA partnership et al.l (120151) gives the 
detailed tuning frequencies for each band. To summarize, the 
1.3 mm and 0.87 mm data (Bands 6 and 7) have four 2 GHz 
(1.875 GHz usable) spectral windows (spws) with 128 chan¬ 
nels per spw (hereafter referred to as a “continuum” spw). 
The 2.9 mm data (Band 3) were taken with two different tun¬ 
ings. One used three continuum spws, along with two narrow 
bandwidth spws centered on HCO + (1-0) and HCN (1-0), each 
with a channel resolution of ~ 0.21 km s -1 and a velocity res¬ 
olution of 0.42 km s -1 . The other 2.9 mm tuning used two 
continuum spws, one narrow bandwidth window on 12 CO (1- 
0), and another spanning four hyperfine transitions of CN (1- 
0). These data have a channel resolution of ~ 0.16 km s -1 
and a velocity resolution of 0.32 km s -1 . The scheduling 
blocks (SBs) were designed to run for 60-70 minutes, with 
approximately 30 minutes on-source per execution, and each 
included about 30 antennas. Each of the 2.9 mm SBs were 
executed seven times, while the 1.3 mm and 0.87 mm SBs 
were executed nine and ten times respectively. The start times 
were staggered to obtain good hour angle coverage for each 
band, and hence good uv-coverage. It cannot be overstated 
how crucial the uv-coverage is to the successful imaging of 
the complex structure of HL Tau. 

The calibration and imaging for the continuum and spectral 
line data are described in § [Al The final synthesized beams 
and rms noise levels for the continuum data are given in Ta¬ 
ble m Cubes with 0.25 km s 1 velocity width channels were 
made for the four observed lines at 1." 1 angular resolution. 
Additionally, an HCO + (1-0) cube with 0. // 25 resolution was 
also created (rms noise levels for the cubes are given in §lAb. 
All of the images described in this paper are publicly available 
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from the ALMA Science Verification paeeFI. 


3.1.2. Continuum Emission 


3. RESULTS 

Figure [It demonstrates the complex outflow emission that 
fills the HL Tau region as traced by both atomic and molec¬ 
ular gas. The strong Ha and [SII] emission (observed with 
the HST) and weak blue-shifted 12 CO (1-0) to the NE of 
HL Tau as well as the strong red-shifted 12 CO (1-0) emis¬ 
sion to the SW (observed with ALMA) are perpendicular to 
the disk, consistent with previous observations (see §1). In¬ 
terestingly, the Ha , [SII], and blue-shifted 12 CO (1-0) emis¬ 
sion to the SE of HL Tau parallel to the disk are difficult to 
reconcile with a simple disk/outflow scenario, suggesting that 
the blue-shifted outflow has broken out of the parental core 
(Monin et al.l lll996l) . or that there is another - as yet unidenti¬ 
fied - driving source. Unfortunately, the 12 CO (1-0) data are 
missing significant flux (due to a lack of short spacings), and 
have insufficient sensitivity in the outer portions of the field 
of view to warrant deeper analysis of its properties. Figs. |T]3, 
and c show zoomed in views of our serendipitous detections of 
XZ Tau (A and B), and LkHa358; no other continuum sources 
above the local 4 a level were detected. 

3.1. HL Tau 

3.1.1. Position and Proper Motion 

The fitted position for HL Tau in each of the ALMA 
images is given in Table |T] The phase calibrator posi¬ 
tions are accurate to < 1 mas and the positions are consis¬ 
tent between the three observed bands to better than 2 mas 
(consistent with dedicated LBC astrometry experiments, see 
lALMA partnership et al.ll2Q15h : thus, we assume 2 mas as the 
ab solute ALMA posi tion uncertainty. The position reported 
bv lKwon et al.1 ( 2011 ) from 1.3 mm CARMA observations is 
04 h 31 m 38 s .418 +18°13 / 57. // 37 (J2000, epoch 2009.08). The 
phase calibrator for CARMA observations (J0510+1800) had 
a position accurate to better than 1 mas, and we assume an 
overall astrometric uncertainty of 5 mas for this measurement. 
Adding the two uncertainties in quadrature, the measured 
angular separation between the CARMA position and the 
1.0 mm ALMA position (epoch 2014.83, TableQ} is AR.A.= 
+106.1 ± 5.6 mas, ADec= -128.0 ± 5.6 mas. With a time 
span of 5.75 years, this separation amounts to a millimeter- 
derived proper motion of vra. = +18.5 =b 1 mas year -1 , and 
VDec. = -22.3 =b 1 mas year -1 . This res ult is in good agreemen t 
with the proper motion adopted by iGuilloteau et al. 12011 ) 
based on 12 years of Plateau de Bure Interferometer (PdBI) 
measurements (+14,-20 mas year -1 ). For comparison, the 
weighted average of optical and infrared proper motions for 
HL Tau listed in VizieS since 2000 (e .g. iDias et all 120141: 
IZacharias et al.ll2013l:lRoeser et al.ll2010l) is v R a. =+2.0 zb 2.4 
mas year -1 , vo e c. = -20.9 zb 2.4 mas year -1 . Apparently, the 
optical/IR proper motion for HL Tau is reasonably accurate in 
Dec., but is significantly underestimated in R.A. This discrep¬ 
ancy is not surprising given that the optical and even near-IR 
emission from HL Tau is dominated by reflection nebulosity 
rather than stellar light (see for example Fig. |T^). 

21 The raw data, calibration scripts (including detailed ab¬ 
solute flux scale information), and images are available from 
http: //a l mascience .org/alma-data/science-verificati 
Additional details along with the full imaging scripts are available from 

http://casaguides.nrao.edu/index.php?title=ALMA2014_ 

22 http://vizier.u-strasbg.fr/viz-bin/VizieR 


Figures |2t,b, and c show the HL Tau continuum at all three 
observed bands. These images reveal for the first time the 
remarkable morphology of the HL Tau protoplanetary disk, 
with a complex pattern of alternating bright and dark rings. 
With the caveats of lower spatial resolution at 2.9 mm and 
reduced phase stability and hence image fidelity at 0.87 mm, 
the agreement between the different wavelengths is excellent. 
Although the antenna configuration was not ideal (baselines 
were concentrated at lengths < 200 m and > 500 m), a rea¬ 
sonably Gaussian synthesized beam was achieved as demon¬ 
strated in Fig. |2]I which shows the 1.3 mm point spread func¬ 
tion (psf). The combined Band 6+7 image at 1.0 mm and its 
corresponding spectral index map are shown in Figs.[2k, and f 
(see §|A]for details). Notably, while the ring structure remains 
clearly visible and much improved over the 0.87 mm data 
alone, the combined 1.0 mm image does not show some of 
the (potentially interesting) azimuthal changes in ring bright¬ 
ness seen in the individual images. These might be due to 
variations in dust opacity, but they could also just be the re¬ 
sult of differences in the uv-coverage at the two wavelengths. 
Thus, we strongly caution against over-interpretation of sub¬ 
tle features unless verified by detailed modeling or future even 
higher fidelity observations. 

The continuum properties of the HL Tau disk are given in 
Tabled] A wide range of flux densities is available in the 
literature. Considering recent interferometric observations, 
2.7 mm flux densities range from 94.1 zb 0.9 to 120 zb 4 mJy 
and 1 .3 mm flux densiti e s range from 700 zb 10 to 818 zb 10.8 
mJy dKwon et all 1201 ll : IGuilloteau et all 1201 ll) . Using the 
more precise 2.7 mm measurement, and scaling by a spectral 
index a = 3, there is good agreement with the 2.9 mm ALMA 
measurement (74.3 zb 0.4 mJy). At 1.3 mm, past observations 
bracket the total flux measured by ALMA (744.1 zb 1.5 mJy). 
At 0.87 mm, only one interferometric measurement is avail¬ 
able: 1300 zb 300 mJy (S MA with 2. // l x U'0 resolution 
iLumbreras & Z apata 20141). but this is 60% smaller than the 
ALMA measurement of 2140.8 zb 3.7 mJy. In contrast, previ¬ 
ous bolometric JCMT measur ements at 0.86 mm wavel ength 
have obtained 2360 zb 90 mJy ([Andrews & William^!20051) . in 
much closer agreement with the ALMA data, especially con¬ 
sidering the difference in angular resolution. This latter com¬ 
parison suggests that only ~ 10% of the emission is resolved 
out by the ALMA data. 

Using the ALMA integrated flux densities from Table [TJ 
the average spectral index is a = 2.77±0.13, with no sig¬ 
nificant curvature evident from 2.9 to 0.87 mm, despite the 
expectation that the emission would be more optically thin 
at the longer wavelength. This result is likely a reflection of 
the fact that the central ~ 10 AU is both significantly brighter 
than the surrounding disk and is possibly optically thick even 
at 2.9 mm. Additionally, there may be weak free-free con- 
tami nation in the central region at 2.9 mm (see for exam¬ 
ple IWilner et~aI1 [19 96: Carrasco-Gonzale z~et all 120091) . For 
comparison,[Andrews & Williams (20051) find a = 2.53 ±0.13 
based on a compilation of primarily bolometer measurements 
between 1.3 and 0.35 mm, suggesting the average spectral in¬ 
dex does eventually flatten at shorter wavelengths. 

The HL Tau disk is viewed on the sky inclined with respect 
to our line of sight (defined as i = 0° for face-on), and rotated 
by an amount defined by the position angle (P.A.; measured 
LBCj^$MJ$orth to East). In order to constrain these parameters we 
first followed the standard practice of fitting a Gaussian to the 
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Fig. 1.— Panel (a) shows an overview of the HL Tau region with red and blue colors mapped to the integrated intensity of the redshifted (7.25 to 20 km s l ) 
and blueshifted (0.0 to 6.0 km s -1 ) 12 CO (1-0), respectively. The ALMA 11 resolution 12 CO (1-0) images (synthesized beam shown in the lower left) have 
not been corrected for primary beam attenuation; the displayed field of view corresponds to the ~ 15% power point. Green shows an HST R-band image (this 
band includes stellar continuum, Ha, and [SII]). The darker diagonal stripe across the image (visible in some viewers) corresponds to a gap in the HST ACS 
WFC detector. Panels (b) and (c) show zoomed in views of XZ Tau (A and B) and LkHa358, respectively. All three panels show ALMA 2.9 mm primary beam 
corrected continuum contours overlaid in blue at 33 x (4, 14,19) n Jy beam -1 (the corresponding synthesized beam is shown in the lower left of panels (b) and (c), 
also see Table[l}- The HST image has been precessed to epoch 2014.83 using the proper motion derived in §3.1.1. 


visibility data. Using baselines < 1000 m 0!'2), we obtain 
i = 46.2 dz 0.2° and P.A.= 138.2 ±0.2° in good agreement with 
past u v-plane estimates dKwon et alJ 120111: iGuillotea u et alJ 
2011). However, the fantastic fidelity and angular resolution 
afforded by the ALMA images allow us to go further, and ex¬ 
plore the properties of the rings independently in the image- 
plane. 

3.1.3. Spatially Resolved Disk Geometry and Spectral Index 

Upon visual inspection, one can plausibly identify seven 
pairs of “bright" and “dark” rings in the 1.0 mm image, we 
label these rings B1..B7 and D1...D7, respectively. As a first 


step, we assumed that the rings trace circular orbits around 
a common center position defined by the 1.0 mm continuum 
peak position (see Table [B, and having the i and PA. derived 
in §3.1.2. Then using a cross-cut along the major axis of the 
image, we determine approximate semi-major axes for the 
fourteen rings. After overlaying these zeroth-order ellipses 
on the image, it was apparent that while approximately cor¬ 
rect, these are not a good fit to the rings in detail. To refine 
the parameters, we defined a discrete set of points along each 
zeroth-order ellipse, at a Nyquist sampled interval with re¬ 
spect to the synthesized beam. Then the position of each point 
was moved to the nearest local radial maximum (or minimum 













ALMA Observations of HL Tau 


5 


TABLE 1 

Measured Continuum Parameters 21 


Source 

Position 15 

rms 

Peak c Flux Density 0 

Size 0 


R.A. (J2000) 

Dec. (J2000) 

{jiiy beam -1 ) 

(mJy beam 1 ) 

(mJy) 

(mas x mas [deg]) 


2.9 mm (Band 3, 101.9 GHz) 85.3 x 61.1 mas,PA=-179° 

HL Tau 

04:31:38.4253 (0.0002) 

18:13:57.240 (0.004) 

9 

4.13(0.03) 

74.3 (0.4) 

- 1750 

XZ Tau A 

04:31:40.09714 (0.00004) 

18:13:56.674 (0.001) 

24 

1.99 (0.04) 

2.71 (0.08) 

<54 

XZ Tau B 

04:31:40.08218 (0.00007) 

18:13:56.845 (0.001) 

24 

1.30 (0.04) 

1.83 (0.08) 

<54 

LkHcx358 

04:31:36.15383 (0.00007) 

18:13:42:919 (0.002) 

26 

1.07 (0.03) 

3.7 (0.1) 

150 x 84 [+170] (6x3 [3]) 

1.3 mm (Band 6, 233.0 GHz) 35.1 x 21.8 mas, PA=+11° 

HL Tau 

04:31:38.42548 (0.00008) 

18:13:57.242 (0.001) 

11 

6.48 (0.03) 

744.1 (1.5) 

- 1960 



0.87 mm (Band 7, 343.5 GHz) 29.9 x 19.0 mas, PA= -176' 



HL Tau 

04:31:38.4254 (0.0001) 

18:13:57.242 (0.001) 

23 

11.56 (0.07) 

2140.8 (3.7) 

- 1960 


1.0 mm (Combined Band 6+7, 287.2 GHz) 33.5 x 21.1 mas, PA: 

= +9° 


HL Tau 

04:31:38.42545 (0.00009) 

18:13:57.242 (0.001) 

12 

9.79 (0.04) 

1441.5(1.8) 

- 1960 


a For each band, the mean continuum frequency and synthesized beam is given after the wavelength. Uncertainties are given in parenthesis after 
each quantity. All measurements are made from images corrected for primary beam attenuation. 

b The peak positions (and uncertainties) were measured using Gaussian fitting, for HL Tau the fit was restricted to the region inside the smallest 
radii gap. 

c For XZ Tau and LkHo;358, the peak intensity, flux density and size (or upper limit) were measured using Gaussian fitting. For HL Tau the peak 
intensity corresponds to the peak pixel value with an uncertainty of 3 a (where a is given in the rms column). The flux density was measured 
within the 4cr contour, and the uncertainty is [No. Independent Beams] 0 5 x 3cr. Uncertainties for these quantities do not include the 5% absolute 
flux uncertainty. The size of HL Tau was estimated from the length of the line passing through the continuum peak at a position angle of +138° 
and where it crosses the 4cr contour on each side of the disk. 


TABLE 2 

Properties of Bright and Dark Rings 21 


Ring 

Name 

Position Offset b 
RA.(mas) Dec.(mas) 

Semi-maior 

(AU) 

Inclination 

(°) 

PA 

(°) 

Dl 

-2.0 ±0.9 

-1.9±0.9 

13.2 ± 0.2 

44.9 ±1.1 

137.2 ±1.8 

Bl 

-0.7 ±0.7 

-4.9 ±0.7 

20.4 ±0.1 

46.0 ±0.5 

138.2 ±0.9 

D2 

-0.9 ±0.5 

-3.0 ±0.6 

32.3 ±0.1 

45.5 ±0.3 

138.1 ±0.5 

B2 

-1.6 ±0.5 

-7.9 ±0.6 

38.1 ±0.1 

46.5 ±0.2 

137.3 ±0.4 

D3 



-42 



B3 



-47 



D4 



-50 



B4 



-55 



D5 

-1.8±0.4 

3.4 ±0.5 

64.2 ±0.1 

45.8±0.1 

138.9 ±0.2 

B5 

-0.8 ±0.5 

7.5 ±0.6 

68.8±0.1 

46.6 ±0.1 

138.8±0.2 

D6 

-4.7 ±0.5 

8.3 ±0.5 

73.7 ±0.1 

47.9 ±0.1 

137.0 ±0.2 

B6 

-8.6 ±0.4 

12.3 ±0.4 

81.3 ±0.1 

46.8±0.1 

137.9±0.1 

D7 



-91 



B7 



-97 




a Rings with a complete set of parameters were calculated assuming circular 
orbits that are tilted by the inclination, and rotated in the sky by the PA. using 
the method described in §3.1.3. Rings with only a Semi-major axis entry were 
estimated directly from the image and have an uncertainty of order 1 AU. 
b Offset positions are with respect to the fitted 1.0 mm peak position in Table[l] 

for dark rings). To avoid regions where the rings become less 
distinct, points were discarded if they moved outside the nom¬ 
inal width of the individual rings (5 to 8 AU). Eight rings re¬ 
tained > 55% of the points, to which we subsequently fit an 
ellipse, including its center position, using _a Markov Chain 
Monte Carlo ( MCMC: [Foreman-Mackev et al.112013l) . The re¬ 
sults are listed in Table [2] with the full range of parameters 
given for the eight most distinct rings, and just the semi-major 
axis for the others. It seems likely that the “gap”, “enhance¬ 
ment”, and “clump” observed in VLA 1.3 and 0.7 c m images 
(ICarrasco-Gonzalez et al.11120091: [Greaves et al.lll2008l) at ~ 10, 
20, and 55 AU along the major axis of the disk correspond 
to the Dl, B1 and the combined emission from the B2 to B4 


rings, respectively. 

The weighted average of the best-fit inclination and position 
angle for the eight fitted rings yields i = 46.72° ±0.05° and 
P.A.= 138.02° ±0.07°, consistent with the constraints found 
for the average disk geometry over large scales. However, the 
best-fit ellipses have their center’s offset with respect to the 
peak of the 1.0 mm emission, as can be seen in the equato¬ 
rial offsets reported in Tabled These offsets are statistically 
significant for all but the innermost ring (Dl). Interestingly, 
the magnitude of the position offset increases with orbital dis¬ 
tance from the center. 

Using the weighted average inclination and P.A., we 
have deprojected the combined 1.0 mm visibility data into 
a circularly-symmetric, face-on equivalent view (see Fig¬ 
ure [3k). We have also extracted cross-cuts at an angle of 
138° from both the 1.0 mm continuum image and the spectral 
index map shown in Figs. [2k, f. These cross-cuts are shown 
in Figs. 03, c. The variation in intensity between the bright 
and dark rings is readily apparent. Considering only the fully 
characterized rings, the largest average intensity contrast is 
between the first pair with Dl being 46% less bright than Bl, 
and the smallest contrast is between the 5th pair with D5 be¬ 
ing 15% less bright than B5. Such a drop in intensity could 
be due to a reduction in dust temperature, column density, or 
grain emissivity, or a combination thereof. This figure also 
demonstrates the general trend of having very high S/N on 
the estimate of a in the (brighter) inner parts of the disk with 
increasing uncertainty as the disk intensity decreases. In¬ 
terestingly, the minimum a does not occur precisely at the 
continuum peak position, but is instead offset by ^ 1.5 AU 
(^ 10.7 mas) along the major axis to the SE. The origin of 
this offset is unknown. 

As shown in Figs. 03, c, each dark ring corresponds to a 
local maximum in the spectral index, while each bright ring 
corresponds to a local minimum. We do not find a gradual 
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Fig. 2.— Panels (a), (b), and (c) show 2.9, 1.3, and 0.87 mm ALMA continuum images of HL Tau. Panel (d) shows the 1.3 mm psf for the same FOV as the 
other panels, as well as an inset with an enlarged view of the inner 300 mas centered on the psf’s peak (the other bands show similar patterns). Panels (e) and 
(f) show the image and spectral index maps resulting from the combination of the 1.3 and 0.87 mm data. The spectral index (a) map has been masked where 
OL/oLerror < 4. The synthesized beams are shown in the lower left of each panel, also see Table [I] The range of the colorbar shown for panel (b), at 1.3 mm, 
corresponds to -2xrms to 0.9x the image peak, using the values in Tablelt] The colorscales for panels (a), (c) and (e) are the same except using the values of 
rms and image peak corresponding to each respective wavelength in Table [J 


decrement of the spectral index with radius reported in other 
lower mass protopl anetary systems, albeit with reduced angu - 
lar resolution (e.g. iGuilloteau et al.ll20TTI: iPerez et al.ll2012l) . 
The central peak, the B1 ring, and the B6 ring, show a spec¬ 
tral index a ~ 2 indicative of optically thick dust emission 
(within 3 <t). The observed 7# (see Fig. [3]l) provides a strict 
lower limit on the physical temperature, at a given angular 
resolution. In the limit that the emission is optically thick, T B 
provides a measure of the physical temperature of the mate¬ 
rial where r ~ 1. The 7# for the optically thick continuum 
peak, and azimuthally averaged values for the optically thick 
B1 and B6 rings are 212.4±0.8 K, 59 =b 3 K, and 24 ± 2 K 
(corresponding intensities are 9.79 ±0.04, 2.48 ±0.12, and 


0.85 ± 0.06 mJy beam -1 ), at radii of ~ 0, 20, and 81 AU, re¬ 
spectively. As shown in Fig. [3ji, the observed radial decrease 
in T b for all the bright rings can be roughly characterized 
by a power law with an exponent of ~ -0.65. For compar¬ 
ison, (1201 j) found T b - 25 K at 55 AU from 

r^j \ ,r resolution PdBI data in good agreement with that pre¬ 
dicted from the ALMA T B analysis. However, the observed T B 
power-law should not be mistaken for a model of the physical 
dust temperature, because much of the disk does not have a 
spectral index consistent with optically thick emission. 

The standard model used for circumstellar disks consists 
of a surface layer of grains directly heated by stellar photons 
which subsequently illuminates and heats the dust in the lower 
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Fig. 3.— Panel (a) shows the deprojected 1.0 mm B6+B7 image of HL Tau 
(see §3.1.3); the angular resolution is 38.6 x 19.3 mas (PA. -20.7°). Rings 
for which a full range of ellipse parameters could be fit are labeled hori¬ 
zontally (solid and dashed lines), while the less distinct rings are labeled 
vertically (dotted lines). Panels (b) and (c) show cross-cuts at PA= 138° 
through the continuum peak of the 1.0 mm continuum and spectral index 
images shown in Fig. l2t,f (positive values of distance correspond to the SE 
portion of Fig. |2^). In panel (c) the grey region delineates the statistical lo 
spectral index uncertainty, it does not account for the absolute flux uncer¬ 
tainty. For the fully fitted rings, panels (b) and (c) show dashed lines for the 
locations of dark rings and solid lines for bright rings. The fitted offsets from 
the continuum peak (Table 2) have been taken into account. Panel (d) shows 
the same cross-cut as panel (b) but on a Planck brightness temperature ( Tb ) 
scale. Panels (b) and (d) are shown on a log scale on the y-axis. The dashed 
curve on panel (d) shows a representative power law for Tb as a function of 
radius with an exponent of -0.65 extending from the B1 peak, it is not a fit. 

layers and the midplane (IChiang & Goldreichll119971) . While 
the surface layer is optically thick to stellar photons, it is op¬ 
tically thin at millimeter wavelengths. Thus, the brightness 
temperature observed by ALMA in the optically thick regions 
should correspond to the physical temperature at the r ^ 1 
surface, which will be somewhere between the surface tem¬ 
perature {T s ) and the midplane temperature (T m ). The ob¬ 
served T b in the central peak and the B1 ring do indeed lie 
between the values of T s and T m at those radii predicted from 
the lKwon et al.ll (120111) “thick-disk" model that best fits the full 
HL Tau SED. Furthermore, the T B of all of the dark rings falls 
below T m , consistent with the dark rings being optically thin. 


The T b in the B6 ring is somewhat below the T m prediction of 
33 K at 81 AU (despite having a ~ 2). One possible inter¬ 
pretation for this behavior is that B6 is not optically thick but 
instead has very low (3. However, we note that the lKwon et all 
( 2011 ) "thin-disk" model, which provides the best match to 
the previous millimeter data alone (due to dust settling to the 
mid-plane), predicts T m values that fall below the observed 
T b (already a lower limit on the physical temperature) in all 
of the dark rings. Clearly, more sophisticated modeling of 
HL Tau’s physical dust temperature distribution that includes 
the radial intensity and spectral index variations, as well as 
radial changes in the scale-height will be fruitful. 

As shown in Fig.[3t,d, the dark rings are not completely de¬ 
void of emission, and the spectral index at their location dif¬ 
fers from the optically thick expectation of a = 2, with larger 
a values ranging from ~ 2.3 to 3.0. This observation rules 
out the possibility that the reduction in emission is due solely 
to radial temperature variations. Instead, these regions are 
likely optically thin, consistent with their lower observed T B . 
An estimate for the value of the dust opacity spectral index (3 
can be de rived in the optically thin l imit for warm dust where 
/ 3 = ol-2 iBeckwith & Sargent! 19911) . We estimate [3 ~ 0.8 for 
Dl, fj ~ 0.7 for D2, /? ~ 0.6 for D5 and /? ~ 0.3 for D6, con¬ 
sistent with some amount of grain growth and evolution inside 
these dark rings. However, further physical modeling of the 
disk dust continuum emission is needed to confirm these es¬ 
timates. Finally, we caution that observations which cannot 
resolve the HL Tau morphology will combine emission from 
optically thick and thin regions, driving any derived values of 
dust opacity spectral index to be lower than reality. 

3.1.4. Compact Spectral Lines 

At 11 resolution, the HCO + (1-0) emission shows a mor¬ 
phology similar to the 12 CO (1-0) shown in Fig. Hk al¬ 
beit somewhat less extended. Interestingly, it also shows a 
barely resolved velocity gradient across the HL Tau disk it¬ 
self, though it is confused with the surrounding outflow gas, 
especially to the NE of the disk. Fortunately, at 0"25 reso¬ 
lution most of the outflow emission is resolved out, making 
it possible to spatially resolve the morphology of the HL Tau 
molecular gas disk for the first time. The 0."25 resolution 
HCO + (1-0) channel maps shown in Figure 0] reveal: (1) A 
notable deficit of blue vs. red-shifted emission; (2) HCO + (1- 
0) absorption near the systemic velocity 6.5-7.0 km s -1 ; (3) A 
roughly Keplerian velocity distribution with a detectable ra¬ 
dial velocity range of 2.0 to 12.0 km s -1 ; (4) Comparable gas 
and continuum disk sizes (at least at the current sensitivity 
level). To further explore points (1) and (2), Figure [5] shows 
spectra from the 11 resolution cubes toward the continuum 
peak. All of the lines show absorption at 6.0 to 7.0 km s -1 , 
indeed CN and HCN (1-0) are only detected in absorption. 
Additionally, the absorption is non-Gaussian in shape, instead 
showing a gradual increase in absorption on the blue-shifted 
side compared to a steeper rise on the red-shifted side. As 
shown in Fig. [1^, the blue-shifted outflow emerges from the 
NE of HL Tau and, based on the continuum disk orientation, 
propagates toward us at i « 47°. Self-absorption by this out¬ 
flowing gas is likely responsible for both the non-Gaussian 
line shape and the deficit of observable blue-shifted disk emis¬ 
sion. The deepest absorption for the CN and HCN transitions 
occurs at a LSRK velocity of 7.0 d= 0.2 km s -1 which we take 
to be the rest velocity of the system. 

Under the assumption of circular Keplerian motion, and 
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Right Ascension (J2000) 

Fig. 4.— Contours show channel maps of the HCO + (1-0) emission with ~ 0725 angular resolution and 0.5 km s -1 velocity resolution overlaid on the 1.0 mm 
B6+B7 continuum image in greyscale. The LSRK velocity of the emission is shown in the upper right corner of each panel. The contour levels are at 1.9(lcr)x 
(-3.5, 3.5, 5.0, 7.0) mJy beam -1 and the color indicates the change from blue- to red-shifted emission. Negative contours are shown dashed. The cross-cut 
position angle of +138° is shown on the 2.5 and 12.0 km s -1 panels and the HCO + (1-0) synthesized beam is shown in the lower left corner of the last panel. 


noting that the velocity extrema of HCO + (1-0) emission oc¬ 
curs at AV =b 5 km s -1 (see Fig. 0]) from the systemic veloc¬ 
ity (7.0 km s -1 ) at a radius of ~ 25 AU, we find that with 
i = 47°, the enclosed mass is ~ 1.3 M 0 . This value is near 
the high end of the range previously reported for HL Tau (see 
§1). However, it is clear that even on the less absorbed red- 
shifted side, the velocity pattern is not so simple, and may 
for example have a c ontribution from infall (see for example 
[Gomez & D’Alessio 2000 ). Future detailed radiative transfer 
analysis coupled with a physical model will be required to 
reproduce the complex HCO + (1-0) absorption and emission 
toward the HL Tau disk in order to obtain a more accurate 
kinematic stellar mass. 

3.2. XZ Tau 2.9 mm Continuum 

At 2.9 mm we resolve the known multiple system XZ Tau 
into two components, A and B (Fig. [Tb, Table \J] also see 
iForgan et~aI1l20l4 ICarrasco-Gonzalez et al.ll2QQ9h separated 
by 273±1 mas at a position angle of 12 8.7±0.5°. This sep- 
aration is 8% smaller than predicted bv IForgan et al . (2014) 
for a circular, face-on orbit (296±1 mas), particularly in R.A. 
This is likely an indication that the orbit is not face-on, but this 
will require future observations to confirm and quantify. Like 
IForgan etal] ( 2014 ). we find no evidence for component “C" 


(a pu tative third star) reported by ICarrasco-Gonzalez et al.1 
<2001 at 7 mm. Using the 2.9 mm flux d ensities from TableQ] 
and the JVLA 7 mm flux densities from IForgan et al.l (120141) . 
we find spectral indices of +1.8 ± 0.5 for both XZ Tau A and 
B, suggesting both have a_free-free co mpon ent in addition to 
dust emission (see also ICarrasco-Gonzalez et al.ll2009l) . 


3.3. LkHa358 2.9 mm Continuum 

At 2.9 mm, we have resolved the LkHa358 disk for the first 
time in millimeter continuum to a size of only 21 ± 1 AU with 
an inclination an_gle of 56 ± 2° at a PA.= 170 d= 3° (see Fig|T]c, 
TablefTT). [Schaefer et al.l(120091) used PdBI to observe 12 CO (2- 
1) with a resolution of 1 ."38 x 0."83 toward this source, ob¬ 
taining an inclination of 28 ± 9°, Vioosin/=1.35±0.04 km s -1 
(the radial velocity at a radius of 100 AU), and a dynami¬ 
cal mass of 0.5-2 M 0 . These authors note that the inclina¬ 
tion angle is not well-constrained by the PdBI data. If we 
assume instead that the molecular gas has the same incli¬ 
nation as the continuum the dynamical mass is significantly 
smaller, 0.3 ±0.1 M 0 . The ALMA 2.9 mm integrated flux 
density measured for LkHa358 (3.7 ±0.1 mJy) is in reason¬ 
able agreement with that measured by ISchaefer et al. ( 2009 ) 
at 2.7 mm: 4.0 ± 0.6 mJy. 
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Fig. 5.— Continuum-subtracted spectra with angular resolution 11 from 
the indicated transitions toward the HL Tau continuum peak. The adopted 
LSRK systemic velocity of 7 km s _l is indicated by the vertical dotted line. 
For comparison the 2.9 mm peak continuum surface brightness at 1 reso¬ 
lution is 5.75 ± 0.01 K. 

4. DISCUSSION 

What is the nature of the remarkable pattern of bright and 
dark rings in the continuum emission from the HL Tau pro¬ 
toplanetary disk? Here we highlight four key observational 
findings. First, the correlated radial structure of higher spec¬ 
tral index and lower brightness temperature in the dark rings, 
relative to the bright rings, demonstrates that the optical depth 
must be lower in the dark rings. Second, the fact that the 
centers of nearly all of the rings are significantly offset from 
the dust continuum peak provides compelling evidence that 
the rings are not circular as assumed, but rather have some 
eccentricity with HL Tau at one focus. Third, the magni¬ 
tude of the offset increases with ring radius, which is tanta- 
lizingly congruent with the observed increase in orbital ec¬ 
centricity for exoplanets at large orbital radii dButler et al.l 
120061: IShen & Turned 120081) . Fourth, several of these rings 
appear to be in resonances with each other. Assuming Kep- 
lerian orbits and neglecting the mass between the rings, the 
first four dark rings have orbital periods in the proportion 
D1:D2:D3:D4 = 1:4:6:8. There also appear to be resonances 
between bright and dark rings, with D2 in a 2:1 resonance 
with B1 and in a 1:4 resonance with B£0. These mean mo¬ 
tion resonances (MMRs) are calculated to be within the con¬ 
fidence intervals obtained on the semi-major axes. Given the 
precise constraints obtained on the orbital radius of the bright 
and dark rings, the n:m resonances found here have uncer¬ 
tainties between 0.01 and 0.09, making these resonance ra¬ 
tios unlikely to arise from random chance. Although multi¬ 
planetary systems show a diverse architectures and many 

23 Other bright rings are close but not exactly at the following integer reso¬ 
nance ratios: D2:B2:B3 are near to a 3:4:5 resonance, while D2:B5 are close 
to a 1:3 resonance. 


pairs of exoplanets are far from being in MMR dLissauer et aTI 
hoTlh . a growing number of stable systems do exist in or near 
MMR (iZhang et al.ll 20141) . Gravitational interactions between 
planets and their parent disk can lead to resonant and near¬ 
resonant systems (see iBaruteau & Papaloizoull2013[ and ref¬ 
erences therein). Collectively, these four independent lines of 
evidence suggest that the observed dark rings are low column 
density gaps in the disk material arising from the process of 
planet formation. 

5. CONCLUSIONS 

In this paper, we present multiwavelength ALMA images 
of the HL Tau protoplanetary disk at resolutions as fine as a 
few AU. We derive precise measurements of the disk inclina¬ 
tion and position angle, as well as the proper motion of the 
central protostar. We also show the first molecular line obser¬ 
vations to spatially resolve the disk kinematics. Under the as¬ 
sumption of Keplerian rotation, the mass of the central star is 
M* ~ 1.3 M©, but further refinement will require detailed ra¬ 
diative transfer calculations. The remarkable pattern of bright 
and dark circumstellar rings in the continuum images exhibits 
a corresponding structure in the spectral index image, reveal¬ 
ing that the dust in the central core and several of the bright 
rings are optically thick, while the dust in the dark rings - 
which are not completely empty - shows evidence of grain 
growth. Several characteristics of these rings, including an 
increase in eccentricity with radius and numerous resonances, 
suggest that the dark rings are gaps arising from the process 
of planet formation. Modeling of the disk temperature and 
density structure in both gas and dust, coupled with planet 
formation processes will be crucial to the goal of reproducing 
HL Tau’s morphology from theory (e.g.. IDong et alJl2014l) . 
These transformational ALMA observations of HL Tau herald 
a new era in the study of protoplanetary disks that promises to 
unearth the architecture of extrasolar multi-planetary systems 
during their epoch of formation, and thereby illuminate the 
origin of our own Solar System. 

This paper makes use of the following ALMA data: 
ADS/JAO.ALMA#2011.0.00015.SV. ALMA is a partner¬ 
ship of ESO (representing its member states), NSF (USA) 
and NINS (Japan), together with NRC (Canada), NSC and 
ASIA A (Taiwan), and KASI (Republic of Korea), in co¬ 
operation with the Republic of Chile. The Joint ALMA 
Observatory is operated by ESO, AUI/NRAO and NAOJ. 
The National Radio Astronomy Observatory is a facility of 
the National Science Foundation operated under coopera¬ 
tive agreement by Associated Universities, Inc. This re¬ 
search made use of Astropy, a community-developed core 
Python package for Astronomy dAstropy Collaboration et al.l 
120131) . as well as the VizieR catalogue access tool, CDS, 
Strasbourg, France.. This paper also makes use of observa¬ 
tions made with the NASA/ESA Hubble Space Telescope, 
and obtained from the Hubble Legacy Archive, which is a 
collaboration between the Space Telescope Science Institute 
(STScI/NASA), the Space Telescope European Coordinating 
Facility (ST-ECF/ESA) and the Canadian Astronomy Data 
Centre (CADC/NRC/CSA). 

Facilities: ALMA. 

24 Note that HR8799 is the only currently known multi-exoplanet system 
spanning a comparable range of _semi-major axes as the gaps in the HL Tau 
disk (see IGozdziewski & Migaszewskil 12014 and references therein). The 
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APPENDIX 

A. IMAGING DETAILS 


The data were calibrated in CASA 4.2.2 using the same procedure adopted for standard science operations, with 
two exceptions. An imprecise position was inadvertently used for the phase calibrator J0431+2037 employed for 
the 2.9 and 1.3 mm data. The shift was small 15 mas), but noticeable when compared to the 0.87 mm data, 
for which a precise position for J0431+1731 was employed. The position of J0431+2037 was corrected to (J2000) 
04 h 31 m 03.76137 s +20°37 / 34.2649 // before calibration, after which excellent position agreement for HL Tau was achieved 
across all three Bands (see fitted positions in Table [TJ. Both final phase calibrator positions were taken from 
http : / /astrogeo . org/vlbi/solutions/rf c_2 012b/rf c_2 012b_cat. txt. The small shifts in position due 
to proper motion over the month long observing span have been ignored. Additionally, the final LBC antenna position correc¬ 
tions were applied before the data were calibrated (see lALMA partnership et all 120151) . For each execution, the bandpass and 
absolute flux calibrators were dynamically chosen by the ALMA online system and alternated between the frequently monitored 
quasars J0423-0120 and J0510+1800, depending on the LST start time 22 . Based on comparison of the resulting calibrated phase 
calibrator flux densities across the many executions for each band, we estimate that the absolute flux calibration is good to within 
5% for the final combined data at each wavelength. 

The continuum images for all three bands were made using multi-frequency synthesis, and the visibility weighting employed 
was mid-way between natural and uniform (i.e., CASA clean parameter robust-O.O). Individually, the 1.3 mm and 0.87 mm 
data have a fractional bandwidth small enough (< 10%) to ignore spectral index effects in the imaging. In contrast, the range of 
continuum spws available for the 2.9 mm image (90.8 to 102.9 GHz), if all combined, would require accounting for at least a linear 
spectral index slope across the band, given the expected range of dust spectral indices a = 2-4 (S u oc i/*). Ho wever , inclusion of 
the spectral slope requires good signal-to-noise ratio (S/N) at both frequency extrema of the calculation dRau & Corn wellll20 111) , 
a condition the 2.9 mm data did not fulfill since four of the five available continuum spws were clustered at ~ 102 GHz. Thus, 
the 2.9 mm continuum image was created excluding the 90.8 GHz spw and without accounting for the spectral index. To make 
an image with the highest possible combination of fidelity and angular resolution we also made a combined Band 6+7 (1.0 mm) 
image (fractional bandwidth 44%) accounting for spatial variations in the spectral index (CASA clean parameter nterms=2). 
In this case, the S/N is high at both extrema and the resulting image and spectral index map are re liable (the derived a have a 
statistical S/N of up to several hundred). Additionally, we used the multiscale imaging option (lCornwellll2008l) for all three bands 
with scales of 0, 5, and 15, approximately corresponding to 0, 1, and 3x the synthesized beam. Using these imaging parameters, 
the combined continuum data at each band was iteratively self-calibrated. The continuum peak position in all bands changed by 
<0.8 mas following self-calibration. The final synthesized beams and rms noise levels are given in Table[T] 

Even with ALMA’s impressive collecting area, the surface brightness sensitivity at very high angular resolution is limited 
for detecting narrow thermal lines. To improve the surface brightness sensitivity, the data for the four Band 3 line transitions 
were tapered to a resolution of ~ 1"1 (and later convolved to exactly this value) and imaged with robust=0.5 , and a velocity 
width of 0.25 km s -1 . The 12 CO (1-0) and HCO + (1-0) lines both show significant extended emission at this resolution and were 
cleaned with the same multiscale parameters as the continuum (an additional scale of 45 was used for CO). The CN (1-0) and 
HCN (1-0) lines are only detected in absorption toward the continuum peak. The resulting rms noise levels for the 11 resolution 
HCN (1-0), HCO + (1-0), CN (1-0), and 12 CO (1-0) cubes are 2.5, 2.5, 5.0, and 9.0 mJy beanT 1 , respectively (the l."l taper 
reduces the effective number of antennas to « 20). Additionally, an HCO + (1-0) line cube was made at a taper corresponding 
to an angular resolution of ~ 0."25 (0."28 x 0."22, PA= +40.5°) with robust=0.0 , and 0.25 km s -1 channels in order to explore 
the compact disk emission; the rms noise in this cube is 2.4 mJy beam -1 . At this resolution 12 CO (1-0) still shows significant 
confusing contamination from outflow emission as well as copious self-absorption, while the CN (1-0) and HCN (1-0) lines are 
only marginally detected. We note that the continuum self-calibration was not applied to the line data, as tests demonstrated that 
this did not improve the S/N. 
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